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ABSTRACT 


An  investigation  is  made  into  the  accuracy  of  the  description  of  the 
actual  mixing-ratio  profile  by  a  power-law  approximation.   The  approxima- 
tion is  obtained  by  a  least-squares  technique  in  order  to  derive  the  best- 
fit  exponent   A  of  the  pressure  profile.   The    A  -values  were  computed 
for  soundings  at  12  geographically  diverse  stations  over  the  United  States 
during  the  period  16-23  March,  1971.   The    A  -values  at  each  station 
were  found  to  undergo  a  time-variation  based  mainly  on  the  synoptic-scale 
variations  in  the  mixing-ratio  profiles.   A  stepwise  multiple  regression 
procedure  involving  up  to  four  variables  from  the  temperature-humidity 
soundings  was  utilized  in  order  to  "predict"  the  value  of   A  from  gross- 
parameters  of  the  soundings  by  station.   Tests  were  performed  at  four 
stations  in  comparing  the  observed  precipitable  water  vapor  W  with  that 
computed  from  the  power-profile  Wcs  which  depended  upon   A  computed  from 
the  sounding.   The  values  of  Wcs  gave  a  high  correlation,  R(WCS,  W  )  ]>  0.97, 
for  each  of  the  four  stations.   However,  the  values  of  W   using   A  from 
the  predictive  multiple  regression  equation  gave  somewhat  larger  standard 
errors  at  each  of  the  four  stations.   The  W   -values  when  compared  to  W 
indicated  some  measure  of  predictive  skill. 
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I.   INTRODUCTION 

Reitan  (1963)  established  a  correlation  factor  of  0.98  between  the 
logarithm  of  the  monthly  precipitable  water  vapor  and  monthly  dewpoints 
based  on  540  observations  at  numerous  locations.   The  regression  formula 
used  by  Reitan  was 

In  W  =  A  +  B  td  (1) 

where  A  =  -0.981  and  B  =  0.0341  are  the  best  fit  coefficients  resulting 
from  the  regression  procedure. 

Smith  (1966)  has  shown  that  an  equation  of  form  (1)  may  be  derived 
from  the  well  known  empirical  equation  of  Tetens  (1930) 

e  =  EQ  x  10<a  fcd  "  B)/(td  +  7  >  (2) 

provided  that  the  mixing-ratio  profile  may  be  parameterized.   The  constants 
are  EQ  =  6108  dynes  cm-2,  a  -  7.5,  B  =  238.1  F,  and   Y  =  395.1  F.   The 
required  parameterization  of  the  mixing-ratio  can  be  accomplished  with 
reasonable  accuracy  by  the  power  law  approximation 

w  =  w0  (p/p0)  (3) 

By  combining  equations  (2)  and  (3)  with 

«o  =  c  C  e0/Po)  (4) 

Smith  was  able  to  show  that 

In  W  »  In  f  EQ/g  +  2.3026  (a  td  -  B)/440  -... 

(5) 
-...-ln(  X  +  1) 

Since  climatological  values  of  the  exponent  A.  are  variable  with  latitude, 

Smith  concluded  that  the  coefficient  A  of  (1)  is  likewise  dependent  upon 

latitude  in  Reitan' s  formulation. 
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TABLE  I 
A.  -values  from  Smith  (1966) 


LATITUDINAL 
ZONE  -°N 

WINTER 

SPRING 

SUMMER 

FALL 

ANNUAL 
AVERAGE 

0  - 

10 

3.37 

2.85 

2.80 

2.64 

2.91 

10  - 

20 

2.99 

3.02 

2.70 

2.93 

2.91 

20  - 

30 

3.60 

3.00 

2.98 

2.93 

3.21 

30  - 

40 

3.04 

3.11 

2.92 

2.94 

3.00 

40  - 

50 

2.70 

2.95 

2.77 

2.71 

2.78 

50  - 

60 

2.52 

3.07 

2.67 

2.93 

2.79 

60  - 

70 

1.76 

2.69 

2.61 

2.61 

2.41 

70  - 

80 

1.60 

1.67 

2.24 

2.63 

2.03 

80  - 

90 

1.11 

1.44 

1.94 

2.02 

1.62 

NORTHERN 

HEMISPHERE 

AVERAGE 

2.52 

2.64 

2.62 

2.70 

2.61 
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Bolsenga  has  attempted  to  adapt  Reitan's  formulation  (1)  to  daily 
and  even  to  hourly  observations  and  finds  that  the  explained  variance  of 
total  water  vapor  falls  off  to  12%   and  64%  in  the  respective  cases.   This 
could  have  been  anticipated  since  total  water  vapor  depends  in  each  case 
upon  the  entire  mixing-ratio  profile,  which  depends  less  upon  surface 
dewpoint  when  a  decreased  time  interval  of  sampling  is  employed. 

The  equivalence  of  equations  (1)  and  (2)  to  the   X  -profile  (3)  holds 
best,  according  to  Smith,  when  applied  to  clima tological  data.   Climato- 
logical  values  of  X  were  calculated  by  Smith  (1966),  from  London's  global 
mean  radiosonde  data  (1957),  with  resultant   X  -values  presented  as  a 
function  of  latitude  and  season  as  shown  in  Table  I. 

The  concept  of  an  accurate  mixing-ratio  moisture  profile  from  sound- 
ings composed  of  few  or  widely  spaced  levels  is  a  very  valuable  tool. 
It  is  especially  important  in  the  field  of  radiation  flux  measurement 
where  values  for  the  mixing-ratio  (moisture  profile  description)  are 
needed  at  all  levels  to  the  top  of  the  atmosphere  but  are  frequently  not 
available  from  sounding  levels  much  above  500  mb.   Recently  Plante  (1973) 
and  Martin  (1973)  have  extended  the  use  of  Smith's  power-law  profile  to 
one  of  the  form 

w  =  w500(p/500)A  (6) 

with  the  reference  or  key  level  at  500  mb  rather  than  a  surface  reference 
level  which  has  a  diurnally  variable  mixing  ratio.   Martin  (1973)  using 
the  Fleet  Numerical  Weather  Central,  Monterey,  California,  computer  print- 
out of  mixing  ratio  at  six  pressure  levels  in  the  vertical  at  FNWC  grid 
points  found  a  correlation  coefficient  of  .998  or  higher  between  observed 
mixing-ratio  and  a  mixing-ratio  computed  by  Eq.  (6),  based  on   X  -values 
calculated  at  each  grid  point  for  the  75  grid  points  considered.   Martin 
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concluded  that  the  power-fit  as  described  by  Eq.  (6)  is  accurate  especially 
with  soundings  having  widely  separated  reporting  levels  in  the  vertical. 
In  this  thesis  the  following  questions  will  be  investigated:   for  a 
given  station  for  a  given  period  of  time,  can  a  power-fit  of  the  nature 
of  (6)  be  obtained  which  is  representative  of  the  actual  moisture  profile 
at  50  mb-level  increments,  and  what  degree  of  accuracy  can  be  expected? 
Further,  the  conditions  which  would  tend  to  invalidate  the  accuracy  of 
the  method  will  be  delineated.   Two  methods  of  verification  of  the  results 
will  be  considered: 

(1)  the  explained  variance  of  mixing-ratio  resulting  from  the  use  of 
power-law  profile,  and 

(2)  the  accuracy  for  the  calculated  precipitable  water  obtained  from 
the  profile-parameterization  compared  to  that  observed. 
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II.   DATA  PROCESSING 

A.   THE  ORIGINAL  DATA 

For  the  water  vapor  analysis,  twelve  stations  were  selected  on  the 
basis  of  their  locations  and  the  availability  of  data.   The  time  chosen 
was  an  eight  day  period  from  16  March  to  24  March,  1971.   The  station 
locations,  elevations,  and  mean  surface  pressures  are  shown  in  Table  II. 
As  can  be  seen  from  Figure  1,  the  stations  selected  lie  in  a  band  across 
the  southern,  and  in  a  second  band  across  the  northern  continental  United 
States.   The  data  was  obtained  from  NWSED,  Ashville,  North  Carolina.   All 
data  levels  were  selected  at  levels  corresponding  to  integral  multiples 
of  50  mb  (up  to  400  mb) ,  in  addition  to  the  surface  level. 

The  following  analysis  was  used  as  the  basis  of  solving  for  X   , 

starting  from  Eq.  (6): 

log  w(p)   =   A  log   p^  (7) 

W500  ^Jtf 

Here  p  ranged  from  the  surface  value  to  p  =  400  mb,  and  corresponding 
values  w(p)  were  computed  at  each  level  p  using  sounding-level  data.   If 
the  left  side  of  Eq.  (7)  is  denoted  as  Y,  then  this  equation  becomes 

Y  -   A  X  (8) 

where  X  =  log  (p/500).   In  the  average  sounding  which  extends  to 
p  =  1000  mb,  there  would  be  14  simultaneous  values  (Y-j_ ,  X-^)  so  that  only 
a  best-fit  solution  for  X   is  feasible. 

This  amounts  to  solving  the  linear  regression  equation,  Eq.  (8),  for 
the  coefficient   X  of  the  independent  variable,  with  the  condition  of  the 
zero-intercept  forced  on  the  constant  term.   In  the  program  library  of  the 
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W.  R.  Church  Computer  Center,  the  most  convenient  program  for  performing 
this  objective  was  BIMED  02R  (Dixon,  1966).   While  this  program  is  primarily 
used  for  stepwise  regression  in  the  case  of  multivariate  regression  analysis, 
it  also  has  the  desirable  option  of  making  the  required  transformations  of 
variables  first  to  ratios,  and  then  to  their  logs  as  in  Eq.  (7).   This 
requires  that  the  sounding  data  has  been  prepared  in  properly  formatted 
form  on  IBM  cards. 

B.   CALCULATION  OF  MIXING-RATIOS  AND  OF  A 

Having  decided  on  the  method  of  handling  the  regression,  the  next 
computation  was  that  of  mixing-ratio  to  be  specified  by  the  power-profile 
fit,  Eq.  (6). 

The  data  as  it  was  received  (by  courtesy  of  Mr.  Russell  Schwanz) ,  was 
card-punched  in  a  form  containing  the  pressure,  temperature  (C),  and 
relative  humidity  at  50  mb-intervals  from  the  surface  to  the  top  of  the 
moisture  sounding  (which  was  generally  at  400  mb) .   The  next  step  then 
was  to  compute  the  mixing-ratio  in  its  approximate  form 


w  =  .62197  (RH)  es  (T)  (9) 

- 

where  es  is  the  saturation  vapor  pressure  computed  over  a  plane  water 
surface  for  1  2  0  C,  and  eg  is  computed  with  respect  to  equilibrium 
over  a  pure,  plane  ice  surface  at  temperatures  T  <  0  C.   In  doing  this  at 
each  sounding  level,  the  two  versions  of  the  Goff-Gratch  equations  for  esw 
and  esi,  respectively,  were  used  in  the  forms  indicated  below  as  Eq.  (10) 
and  Eq.  (11). 

logio  esw  -  -7.90298  (Ts  -1)  +  5.02808  logio  (Ts)-  •••  (10) 

T~  T~ 

...  -1.3816  x  10-7  (ioll.334(l-T/Ts)  -1)  +  . . . 
...  +8.1328  x  (10-3.49149 (Ts/T-l)-i)  +  logio  e^s 
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loglOesi  =   -9.09718(To/T-l)-3.56654  log10(To/T)  +  ...  (11) 

...  +0.876793(l-T/To)  +  log10eio 
where  Ts  =  373.16  K  and  TQ  =  273.16  K. 

In  (10)  esw  is  the  saturation  vapor-pressure  of  pure  water  at  the  steam- 
point  temperature  and  esi  is  that  of  pure  water  at  the  frost-point 
temperature  (T  <  0  C) .   Then  w(p)  is  evaluated  from  either  equation  (10) 
or  (11)  utilizing  Eq.  (9).   Once  the  values  for  w  were  obtained  at  each 
level  in  each  sounding  the  best-fit  value  for  A   was  obtained  in  the 
manner  suggested  in  connection  with  Eq.  (8)  and  Section  III  (A) .   Thus  a 
unique  value  of   A  resulted  for  each  sounding.   The    \  -values  for  each 
station  were  then  grouped  with  other  synoptic  sounding  data  for  each 
individual  sounding  in  order  to  account  for  the  reduced  variance  of  A 
for  the  sounding-sample  for  the  particular  station.   In  this  connection, 
the  stepwise  regression  program  BIMED  02R  was  used  to  determine  the 
explained  variance  of   A  at  each  station.   Comparisons  were  made  of 
computed  w-values  resulting  from  Eq.  (6)  versus  observed  w-values  over 
the  sounding  station. 

C.   CALCULATION  OF  PRECIPITABLE  WATER  VAPOR 

Precipitable  water  was  obtained  by  integration  for  both  the  observed 
soundings  and  for  the  corresponding   A  -profile  cases  and  a  comparison 
was  made  for  several  stations.   The  formulas  used  in  the  calculations  are 
shown  as  follows. 

a.   The  observed  value  of  precipitable  water  vapor  is 


W  - 
P 


P14 

/»  dp  =  _1_  F"y*wi  +  w2  A   P  +  w2  +  w^  A    P  +.  .  .  (12) 

g      980  \J~*         2     1.5        2      2.5 
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.   +  w12  +  w13      A  P  +     w-L3  +  w-^  A 

2  12.5  2 


13.5 


W   = 


b.      The  model   value   of   precipitable  water  vapor   is 
Ps  P 

/w500      ( -2-\  1  dp   =  w5QQ    (5Q0) 

\  500/  g  g    (A+   1) 

ptop 


V500/ 


400 


-  "500 
980 


(500)    |7_is\  A  +1-/400V+1 
(  A   +  1)1  \500/  \500/ 


(13) 


The  comparison  of  the  sample  sets  of  the  two  precipitable  water  values 
will  be  used  as  a  second  measure  to  determine  the  accuracy  of  the   A.  - 
profile  method. 
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III.   STATISTICAL  CONCEPTS 

A.   THE  LEAST  SQUARES  FIT 

In  this  thesis,  the  least-squares  fit  was  employed  in  at  least  two 
contexts.   First,  it  was  employed  in  solving  Eq.  (8)  for   A  at  each 
station  and  sounding  time.   Here  the  variables  Y  and  X  are  given  point- 
wise  by 

X.  =  log   P   and  Y.  -  log  w.  (14),  (15) 

l     °  l     ° 

500  — - — 

"500 
with  i  representing  a  level  numbered  from  one  at  the  surface,  to  two  at 
the  next  higher  level  (where  P  first  becomes  an  integral  multiple  of 
50  mb) ,  etc.   Standard  texts  in  statistics  (e.  g.  Brownlee,  1960)  show 
that  the  parameter  A   of  the  power-profile  as  in  Eq.  (8)  can  be  deduced 
from  the  best-fit  condition 

<£X±2  (16) 

If  the  zero-intercept . had  not  been  required  in  (16),   A  would  have  been 

estimated  from 

A  =  1  S(Xt  -  X)(Yj  -  Y) 

N         o~^2  (17) 

and  the  constant  term  of  the  best-fit  simple  regression  would  then  be 

N  (18) 

where  N  is  the  number  of  sounding-level  observations  entering  into  the 
power-profile  (normally  here  N  =  14) . 
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The  combination  of  conditions  (17)  and  (18)  gives  the  standard  best- 
fit  solution  for  the  simple  independent  variable  case 

Y  =  A  +   AX  (19) 

Equation  (16)  gives  the  corresponding  result  when  the  zero-intercept 
(Yj_  =  X^  =  0)  is  required.   In  either  case  \      is  related  to  the  simple 
correlation  coefficient  R^2  an^  tne  variables  Y  and  X  through  the  standard 
formula 

&12  =   A  _^*  (20) 

Hence  a  value  of  Ry|x  close  to  unity  means  an  effective  least-squares  fit 
by  the  power-profile  and  simultaneously  a  reliable  value  of  A   .   In 
Eq.  (20),  the  statistical  parameters   ff„  and  o    „   are  defined  as 

2  2  2 


a. 


"y 


Yj,  +  y2  +  ...  YNy- 


,       ox=    /xjUxjU  _  xgV 


(21) 


Note  that  ay   an^  °\   °f  (21)  are  applicable  to  the  zero  -intercept  case 
and  apply  to  the  more  general  case  of  the  standard  regression  formula 

Yi  -  Y  =  Ri2_^Y  /x.  -  Xj  (22) 

when  Y-l  is  replaced  by  Y^  -  Y  and  X^  is  replaced  by  X^  -  X. 

It  was  possible  to  apply  BIMED  02R  regression  procedure  to  the 
14-values  of  X^  and  Y^  per  sounding  (that  is,  with  each  sounding  listed 
on  14  cards)  and  to  generate  not  only  the   A.  -value  for  that  sounding, 
but  also  the  simple  correlation  Ri2>  which  is  a  measure  of  the  "goodness 
of  fit"  of  the  power-profile  for  the  particular  sounding  case.   Such 
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results  for  A.   and  R^2  nave  been  tabulated  under  each  station  and  for 
each  observation  time  (Tables  III  through  XIV). 

B.   THE  STEPWISE  REGRESSION  PROCEDURE 

If  there  are  N  samples  of  the  statistically  related  variables 

(Yli>x2i>X3i>  '  *  *  'Xki)  *■  =  !>-'->N'  *-fc  i-s   possible  to  apply  the  least- 
squares  concept  to  minimize  the  least-squares  form 

k 

X^  (Yli  -  Al   -  A2X2i  "•  • -"AkXki)2=  minimum    (23) 

i=l 

The  procedure  of  selecting  the  least-squares  fit  of  coefficients  to 
minimize  Eq.  (23)  is  called  multiple  linear  regression.   In  this  process, 
there  is  no  assurance  that  the  variables  X2,....,X,  used  to  explain  the 
variation  of  Yi  are  the  most  effective,  and,  in  fact,  some  of  the  X^  may 
not  contribute  to  the  explained  variance  in  any  significant  manner. 
BIMED  02R  is  a  statistical  procedure  for  developing  a  multiple 
regression  equation  in  which  the  variable  X.  added  at  each  step  accomplishes 
the  greatest  reduction  in  the  unexplained  sum-of -squares ,  y     (Y.-Y)  .   If 
k-1  variables  are  to  be  added  in  a  stepwise  manner  the  first  one  added 
is  that  which  has  the  largest  simple  correlation  coefficient  R-j^  of  all 
k-1  possible  comparisons.   As  a  result,  the  first  estimate  to  the  desired 
multiple  regression  equation  may  be  obtained  in  the  form: 

Yi  =  aU  +  aL2X2  (24) 

with  a^,  a-^2  selected  as  in  (17),  (18).   A  property  of  this  first-step 
selection  is  that  the  sum-of -squares  explained  by  Eq.  (24)  is  given  by 

(S.S.  Expl)i  =  R12    (Yn-Y]j  +  ...  +  (YiN  -Yi)2     (25) 

A  result  which  is  proved  in  standard  statistics  texts  (e.g.   Crow  et  al, 
1960). 

22 


The  BIMED  02R  program  then  forms  the  data-sample  residuals.  After  one 
step,  all  values  have  the  form 

Yli-ail-a12X2i  =  YU  (!>  (26) 

(1) 
The  program  then  examines  all  of  the  simple  correlations  of  Yi   on  X3,..., 

Xk  (the  remaining  independent  variables)  and  selects  that  which  has  the 

(1) 
highest  simple  correlation  coefficient  among  R(Y^   >X^),  i  =  3,...,k. 

Since  this  simple  correlation  is  based  upon  the  residual  after  the  effect 
of  variable  X2  has  been  removed  (see  Eq.  (26)  ),  it  is  typically  called 
the  first  order  partial  correlation  with  the  notation  Ri3.2>  tne  under- 
standing being  that  the  variable  having  this  highest  correlation  at 

step  2  and  thus  being  selected  is  called  Xo  .   Then  the  regression  equation 

(1) 
relating  Yi   to  X3  is  obtained  as  described  in  connection  with  Eqs.  (21), 

(22),  (23)  and  assumes  the  form 

(1) 

Yl  =  bll  +  b13  X3  (27) 

As  a  result  of  the  second-step  simple  regression,  Eq.  (27),  the  sum-of- 
squares  of  Yi  (not  Yi   )  is  now  statistically  explained  in  the  amount 


(SS  Expl)  1>2  =  Ri>23   (Yn  -Yx) 


_2  _   2 

+  ...  +  0?1N-Yl) 


'N  -  1     (28) 


2 
The  statistical  parameter  R-^  2o  is  the  so  called  multiple  correlation  at 

step  2,  that  is  to  say  of  variableYi  depending  /jointly  on  the  variations 
of  both  X2  and  X3.   These  independent  variables  were  selected  according 
to  the  criteria  of  having  the  largest  simple  correlation,  and  the  largest 
first-order  partial  correlations  respectively. 

The  multiple  correlation  coefficient  satisfies  the  important  chain- 
rule  of  multiplication 

1  -  r123   =    (1-R^2)  (I-R13.2)  <29) 
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and  after  k-1  steps,  this  rule  assumes  the  form 

1  "  R1.23-'-k  "  CI  -  R12)  (1  "  R13.2---k)"-(1  "  Rlk.23'--k-l)   <3°) 

2 
Just  as  Ri  23  describes  the  fractional  explained  sum-of -squares  of  Y  in 

2 
Eq.  (28),  so  1-R-L  23  describes  the  fractional  unexplained  sum-of -squares 

about  the  regression  surface.   This  parameter  has  decreased  relative  to 

2        2 
that  of  step  one,  provided  R-^  23  >  ^12"   ^he  standard  error  of  estimate 

is  related  to  the  multiple  correlation  coefficient  by 

2         2      2 
SY|X  =    °  Y  (!-R12)  (31) 

The  parameter  standard  error,  Syi^j  will  be  listed  in  Tables  III,..., XIV. 
In  (31),  the  subscript  YIX  denotes  the  best  fit  of  Y  on  all  X^. 

The  most  attractive  feature  of  the  stepwise  regression  is  embodied  in 
the  chain-rule  of  Eqs.  (29),  (30),  with  each  step  making  the  selection 
of  the  highest  partial  correlation  coefficient  for  use  in  Eq.  (30).   The 
result  is  that  when  the  unexplained  fractional  sum-of -squares  does  not 
become  significantly  smaller,  the  process  of  adding  variables  to  the  right- 
hand  side  of 

Yl  =  Al  +  A2X2  +  A3  X3  +  • ' •  AkXk  (32) 

should  be  terminated  at  the  (k-l)th  step. 

In  Section  IV,  application  of  BIMED  02R  will  be  used  on  the  16-sample 
values  of  A   derived  at  each  individual  station.   This  is  done  since  it 
turns  out  that  a  time  variation  in  A   occurs  with  the  synoptic  situation 
at  individual  stations.   In  some  cases  this  variation,  as  indicated  by 
the  standard  deviation  O^    ,  is  so  small  as  to  be  limited  by  the  accuracy 
of  the  data.   In  other  cases,  a  rather  large  variation  of  A  occurs  and 
it  was  desired  to  consider  how  the  time-varying  structure  of  the  soundings 
contribute  to  the   A  -variation.   For  this  further  study  each   A  -value 
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was  added  on  an  IBM  card  containing  certain  relevent  parameters  of  the 
sounding  from  which  X  was  derived.   For  all  soundings  and  stations  a 
set  of  six  independent  parameters  (for  example  T(850)  and  w(850)  and 
others  to  be  described  in  Section  IV)  were  included  along  with  the   X  - 
value  for  the  identical  sounding. 

Thus  BIMED  02R  was  also  used  in  the  important  stepwise  sense  in 
describing  the   X  -variation  in  the  sounding  structure  over  the  set  of 
the  16  soundings  per  station.   The  detailed  discussion  of  the  variables 
used  in  the  regression  for  A   is  given  in  Section  IV  (B,  C),  together 
with  the  discussion  of  the  results. 

C.   ESTIMATES  OF  STATISTICAL  SIGNIFICANCE 

1 .   Of    X  -Profiles  at  Sounding  Stations 

Here  the  correlation  coefficient  R-,~  of  the  regression  equation 
(7)  is  the  basis  of  the  estimation  of  the  accuracy  of  the  X  -profile. 
The  significance  of  the  test  is  based  upon  Eq.  (25)  which  gives  the 

result 

2 
Fractional  Explained  Variance  =  L^,  (33) 

whereas  the  fraction  residual  (or  unexplained)  variance  is  given  by 

2 
Fractional  Unexplained  Variance  =  1-R-i?  (34) 

The  significance  test  is  based  upon  the  F-statistic  (which  is  listed  in 

the  BIMED  02R  output)  2 

F(l,  13)        12 (35) 

(1-R2  )/13 
12 

The  critical  value  of  F(l,13)  at  the  99%  confidence  level  is  listed  in 
statistical  tables  as 
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F  01(1,13)  -  9.07 
This  99%  confidence  level  is  satisfied  whenever  R^2  of  Eq.  (34)  exceeds 

R12  (critical  )  ^  .641 
In  all  but  one  of  the   X  -profiles,  R-^  exceeded  this  value,  and  in  this 

exception  R^2  was  0.614,  corresponding  to  a  confidence  level  of  97.5%. 

2 
In  summary,  the  statistic  R-^2  listed  for  each  sounding  at  the  12 

stations,  is  the  explained  fraction  of  the  variance  of  Y  in  that  sounding 

(see  Tables  III,  ...,  XIV). 

Finally,  BIMED  02R  lists  also  the  standard  error  of  the  coefficient  X 

Standard  Error  of  A   =  A  -A/  1-R12     ;  N  =  14  (36) 

This  statistic  appears  in  column  3  of  Tables  III,...  XIV,  and  describes 
the  sharpness  of  the    \  -fit  in  Eq.  (7). 

2.   Significance  Tests  of  the  Temporal  Variation  of   A 

Here  the  fractional  explained  variance  may  be  computed  from 

Eq.  (28)  after  the  kth  entry  of  an  independent  variable  X  drawn  from  the 

k 

synoptic  sounding  variables  has  been  included  in 

A  -  A0+  ALXX  +  ...  +  AkXk  (37) 

The  fractional  explained  variance  is  given  by  R^  in  accordance  with 
(28),  whereas  the  fractional  residual  variance  is  1-R  .   In  general  when 
0"  (^)  is  sizeable,  it  is  possible  to  generate  a  multiple  regression  equation 
of  form  (37),  which  significantly  reduces  the  residual  variance.   In 
general  the  criterion  used  for  entering  additional  variables  for  explaining 
the   X  -variance  was  based  on  a  stepwise  selection  which  added  at  least 
2%  to  the  total  explained  fractional  variance.   This  strategy  was  used 
because  the  F-statistic  for  this  case  is 
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2       2 

Fk  (1,  N-k-1)  =  ^ ± (38) 

(1  -  RR  )/(N  -  k  -  1) 

and  here  the  critical  value  of  F,  is  difficult  to  surpass  with  N  =  16 
observations.   The  added  fractional  explained  variance  at  step  k  is 
identical  to  the  numerator  of  Eq.  (38),  and  is  required  to  be  at  least 
2%  at  each  step. 
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IV.   RESULTS 

A.   DISCUSSION  OF   A  -PROFILES 

This  particular  phase  of  the  research  began  with  computer  cards 
containing  the  following  individual  sounding  data:   pressure  (mb) ,  tempera- 
ture (K) ,  and  mixing  ratio  (gm/kg) .   There  were  16  soundings  of  up  to  14 
levels  at  each  of  the  12  stations  (Table  II). 

The  BIMED  02R  program  was  utilized  in  a  least-squares  context  to 
solve  Eqs.  7  and  16  for  a  best  fit  A.  for  each  sounding  at  each  station. 
The  "zero-intercept"  was  chosen  in  the  correlation  so  that  the  computer 
solution  to  Eq.  (7)  was  of  the  form  of  Eq.  (8)  instead  of  Eq.  (19).   A 
discussion  of  the  statistical  methods  and  parameters  is  given  in 
Section  III  (A).   Table  III  lists  the  best-fit   A.  -value  for  each  avail- 
able sounding  at  Waycross,  Ga . ,  over  the  test  period. 

It  should  be  noted  that  the  least-square  procedure  for  determining 
from 

Y.  =  A  X.  +  e   .  (39) 

ill 

2 
minimizes  the  sum-of -squares  of  the  residuals   £  e   ..   Hence  the  estimator 

Y  may  be  defined  as 

Y  =   A  log  P/500  =  log  wc(P)  (40) 

w500 
where  the  observed  w^qq  appears  in  the  denominator  of  Eq.  (40)  by  the  zero- 
intercept  concept.   On  the  other  hand,  the  observed  Y  is  defined  as 

Y  =  log  w/w^qq  so  that  the  residual  becomes 

et  =  Y  -  Y  =  log  wc(P)  (41) 

w(P) 


28 


In  Table  III,  successive  values  of   A  obtained  by   A  -profile  fit 
[Eq.  (39)1  have  been  listed  along  with  the  standard  error  of  each   A 
determination.   The  small  standard  error  of  A   indicates  the  close  fit 

of  the  A   values  by  Eq.  (16). 

2 
In  addition,  the  fractional  explained  variance  Rio  for  each  sounding 

is  listed  alongside  R-^  under  column  4.   Finally,  the  standard  error  of 
estimate  is  listed  in  the  final  tabular  column  as  defined  by  Eq.  (31). 
The  value  of  Syix  is  related  to  the  standard  deviation  tfy  through 

SYIX  =  C1  "  R12)^*Y  <4?) 

These  ratios  of  Eq.  (42)  obtained  by  applying  the  values  of  Sy|X  and 
from  columns  5  and  6  of  Tables  III  through  XIV  are  generally  of  the 
order  of  0.3  or  lower.   There  are  exceptions,  but  these  can  be  identified 
by  the  values  of  Rio  (generally  Rio  ^  .95)  for  the  respective  soundings. 
In  Table  III,  the  standard  error  may  be  written  in  terms  of  the  mean 
residuals  -  squared  as 

SY|X      £   |logZ£l    /  log  w  +   [Awcl      (43) 


[2   /log  wcV   I  =    log  w 

i=1   \    —)±     I   13  J 


w 


where  the  final  term  in  Eq.  (43)  represents  the  mean  scatter  about  the 

A 

graph  of  Y  =   XX. 

For  example,  at  0000  GMT,  21  March,  1971,  the  standard  error  at  Way- 
cross,  Ga.,  was  0.1131.   Therefore,  .1131  =  log  w  +   |AWCI  ,  hence  it 

w 
follows  that  w  +  |  Awc 1  ■  1.298.   This  means  that  the  mean  residual  ratio, 

w 
after  regression,  for  this  case  is    |Awc'   ~  30%.   Analysis  of  the  actual 

residuals  for  the  soundings  shown  indicates  that  the  majority  of  this  error- 
ratio  occurs  at  the  highest  altitudes  recorded  in  the  sounding  where  w  is 
the  smallest  and  where  actual  instrument  accuracy  is  doubtful. 
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TABLE    III 

Values   of      A     and    related    statistics 
for    the   data-sample    for  Waycross,    Ga . 


DATE 
DAY -GMT 

A -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

A -VALUE 

CORRELATION 
COEFFICIENT 

Rl2   Ri2 

STANDARD 
DEVIATION 

try 

STANDARD 
ERROR 

SYIX 

16-0000 

2.9414 

0.2585 

.953 

.908 

0.5861 

0.1837 

16-1200* 

1.5298 

0.4277 

.704 

.496 

0.4130 

0.3043 

17-0000 

4.5798 

0.1950 

.988 

.976 

0.8809 

0.1387 

17-1200 

2.7855 

0.1246 

.987 

.974 

0.5373 

0.0888 

18-0000 

2.2575 

0.1187 

.983 

.966 

0.4377 

0.0846 

18-1200 

2.6057 

0.1371 

.983 

.966 

0.5057 

0.0978 

19-0000 

2.3259 

0.3495 

.879 

.773 

0.5038 

0.2491 

19-1200 

1.8176 

0.2879 

.868 

.755 

0.3981 

0.2049 

20-0000 

3.6279 

0.2036 

.980 

.960 

0.7030 

0.1447 

20-1200 

2.5391 

0.2308 

.950 

.903 

0.5083 

0.1643 

21-0000 

3.0887 

0.1589 

.983 

.966 

0.5976 

0.1131 

21-1200 

2.4250 

0.1087 

.987 

.974 

0.4680 

0.0775 

22-0000* 

0.9815 

0.3252 

.642 

.412 

0.2911 

0.2316 

22-1200 

3.2945 

0.2850 

.955 

.912 

0.6572 

0.2031 

23-0000 

4.3912 

0.2107 

.985 

.970 

0.8473 

0.1499 

23-1200 

4.5008 

0.2425 

.982 

.964 

0.8713 

0.1724 

TIME  MEAN 

2.8557 

.869 

STD.  DEV. 

1.0401 

*t<  *  > 

Corresponds    to   R12    <   .85. 
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TABLE  IV 

Values  of   A  and  related  statistics 
for  the  data-sample  for  Charleston,  S.  C 


DATE 
DAY-GMT 

A -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

A  -VALUE 

CORRELATION 
COEFFICIENT 

Rl2    Rl2 

STANDARD 
DEVIATION 

CTY 

STANDARD 
ERROR 

SY|X 

16-0000 

2.3199 

0.1174 

.984 

.968 

0.4485 

0.0863 

16-1200 

4.33721 

0.2933 

.972 

.945 

0.8494 

0.2088 

17-0000 

3.9920 

0.1322 

.993 

.986 

0.7646 

0.0941 

17-1200 

3.8305 

0.1629 

.989 

.978 

0.7379 

0.1160 

18-0000 

2.4849 

0.1724 

.970 

.941 

0.4883 

0.1229 

18-1200 

2.7495 

0.1236 

.987 

.974 

0.5319 

0.0883 

19-0000 

2.9405 

0.1308 

.988 

.976 

0.5684 

0.0934 

19-1200* 

1.5187 

0.4050 

.721 

.520 

0.4013 

0.2886 

20-0000 

4.5054 

0.2364 

.983 

.966 

0.8711 

0.1680 

20-1200 

2.7135 

0.2153 

.961 

.924 

0.5369 

0.1532 

21-0000 

4.7245 

0.1834 

.990 

.980 

0.9077 

0.1305 

21-1200 

2.4194 

0.1502 

.976 

.953 

0.4727 

0.1071 

22-0000 

MISSING  : 

DATA 

22-1200 

2.5801 

0.2807 

.931 

.867 

0.5281 

0.2001 

23-0000 

4.4411 

0.2290 

.983 

.966 

0.8597 

0.1631 

23-1200 

3.5360 

0.1560 

.988 

.976 

0.6809 

0.1110 

TIME  MEAN 

3.2729 

.928 

STD.  DEV. 

0.9832 

V  A  ) 

^Corresponds  to  R-^  <  -85 


31 


TABLE  V 

Values  of  A   and  related  statistics 
for  the  data-sample  for  Lake  Charles,  La. 


DATE 
DAY -GMT 

A.  -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

X -VALUE 

CORRELATION 
COEFFICIENT 
2 
R12    R12 

STANDARD 
DEVIATION 

CTY 

STANDARD 
ERROR 

SYIX 

16-0000 

2.9224 

0.1461 

.984 

.968 

0.5653 

0.1041 

16-1200 

4.0954 

0.1065 

.996 

.992 

0.7837 

0.0759 

17-0000 

4.5898 

0.1810 

.990 

.980 

0.8835 

0.1291 

17-1200 

3.5190 

0.2167 

.976 

.953 

0.6878 

0.1547 

18-0000 

4.0538 

0.1916 

.986 

.972 

0.7843 

0.1367 

18-1200 

3.8884 

0.1964 

.984 

.968 

0.7528 

0.1399 

19-0000 

2.7501 

0.0960 

.992 

.984 

0.5273 

0.0684 

19-1200 

3.0506 

0.2307 

.965 

.931 

0.6030 

0.1646 

20-0000 

2.4835 

0.0959 

.990 

.980 

0.4782 

0.0684 

20-1200 

2.8351 

0.1944 

.971 

.943 

0.5575 

0.1388 

21-0000* 

1.1751 

0.2333 

.813 

.661 

0.2757 

0.1665 

21-1200 

1.6059 

0.2349 

.885 

.783 

0.3462 

0.1676 

22-0000 

4.3386 

0.1766 

.992 

.984 

0.7777 

0.1041 

22-1200 

5.2771 

0.1789 

.994 

.988 

0.8705 

0.1016 

23-0000 

4.4083 

0.2264 

.983 

.966 

0.8532 

0.1612 

23-1200 

4.5449 

0.3695 

.963 

.927 

0.8417 

0.2375 

TIME  MEAN 

3.4711 

.936 

STD.  DEV. 

1.1352 

"t(  *  ) 

^Corresponds  to  R12  <  -85. 
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TABLE  VI 

Values  of  A.  and  related  statistics 
for  the  data-sample  for  El  Paso,  Texas 


STANDARD 

CORRELATION 

STANDARD 

STANDARD 

ERROR  OF 

COEFFICIENT 

DEVIATION 

ERROR 

DATE 

A.  -VALUE 

R12 

2 
R12 

SYIX 

DAY -GMT 

CALCULATED 

A  -VALUE 

aY 

16-0000 

3.5147 

0.3078 

.964 

.929 

0.5391 

0.1509 

16-1200 

3.3290 

0.2334 

.976 

.953 

0.5050 

0.1147 

17-0000 

3.4696 

0.1695 

.988 

.976 

0.5195 

0.0832 

17-1200 

3.2455 

0.3847 

.936 

.876 

0.5123 

0.1886 

18-0000 

3.7536 

0.2578 

.997 

.994 

0.5658 

0.1260 

18-1200 

3.3012 

0.1721 

.987 

.970 

0.4938 

0.0843 

19-0000 

3.3597 

0.1046 

.995 

.990 

0.4990 

0.0513 

19-1200 

2.4246 

0.2875 

.936 

.876 

0.3841 

0.1414 

20-0000 

2.5997 

0.2698 

.950 

.902 

0.4057 

0.1327 

20-1200- 

0.6696 

0.1802 

.778 

.605 

0.1311 

0.0868 

21-0000 

2.6959 

0.3116 

.939 

.882 

0.4232 

0.1524 

21-1200 

3.1878 

0.1860 

.983 

.966 

0.4787 

0.0911 

22-0000 

3.3213 

0.2029 

.982 

.964 

0.4995 

0.0994 

22-1200 

2.4686 

0.3511 

.912 

.832 

0.4003 

0.1722 

23-0000 

2.3344 

0.1755 

.973 

.947 

0.3545 

0.0860 

23-1200* 

0.7583 

0.4089 

.526 

.277 

0.2192 

0.1965 

TIME  MEAN 

2.7771 

.871 

STD.  DEV. 

0.9147 

M  A  ) 

"^Corresponds  to  R^2  <  »85 
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TABLE  VII 

Values  of   A.  and  related  statistics 
for  the  data-sample  for  San  Diego,  Calif. 


DATE 
DAY -GMT 

A  -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

X -VALUE 

CORRELATION 
COEFFICIENT 
2 

R12    R12 

STANDARD 
DEVIATION 

STANDARD 
ERROR 

SYIX 

16-0000 

3.5187 

0.0897 

.  996 

.992 

0.6713 

0.0637 

16-1200 

3.8860 

0 . 1442 

.991 

.982 

0.7444 

0.1024 

17-0000 

3.4676 

0.0874 

.996 

.992 

0.6606 

0.0620 

17-1200 

3.2899 

0.1467 

.988 

.976 

0.5932 

0.0943 

18-0000 

3.9886 

0.1914 

.987 

.974 

0.7205 

0.1230 

18-1200 

4.0171 

0.1426 

.992 

.984 

0.7684 

0.1012 

19-0000 

3.6748 

0.1190 

.993 

.986 

0.7023 

0.0845 

19-1200 

2.4739 

0.1570 

.975 

.951 

0.4817 

0.1115 

20-0000* 

1.9672 

0.3994 

.807 

.651 

0.4627 

0.2836 

20-1200 

3.9124 

0.1845 

.986 

.972 

0.7529 

0.1309 

21-0000 

3.3273 

0.2063 

.976 

.953 

0.6470 

0.1465 

21-1200 

3.8186 

0.1946 

.984 

.968 

0.7370 

0.1382 

22-0000 

2.6873 

0.2288 

.956 

.914 

0.5338 

0.1626 

22-1200 

2.3642 

0.2423 

.938 

.880 

0.4786 

0.1721 

23-0000 

2.1203 

0.2521 

.919 

.845 

0.4378 

0.1790 

23-1200 

2.8397 

0.1020 

.992 

.984 

0.5434 

0.0725 

TIME  MEAN 

3.2096 

.938 

STD.  DEV. 

0.7012 

«rt(  X  ) 

^Corresponds  to  R^2  <  .85 
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TABLE   VIII 

Values   of     X     and   related    statistics 
for    the   data-sample    for   Caribou,    Maine 


DATE 
DAY-GMT 

A  -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

A -VALUE 

CORRELATION 
COEFFICIENT 
2 
R12   R12 

STANDARD 
DEVIATION 

STANDARD 
ERROR 

SY|X 

16-0000* 

1.3979 

0.2714 

.841 

.707 

0.2937 

0.1661 

16-1200 

4.3499 

0.0839 

.998 

.996 

0.7705 

0.0534 

17-0000 

3.5403 

0.1644 

.988 

.976 

0.6588 

0.1047 

17-1200 

2.9646 

0.1566 

.984 

.968 

0.5350 

0.1003 

18-0000 

3.1736 

0.2673 

.960 

.922 

0.5869 

0.1711 

18-1200* 

1.5775 

0.4351 

.723 

.523 

0.3878 

0.2788 

19-0000 

4.6004 

0.1396 

.995 

.990 

0.8453 

0.0884 

19-1200 

3.7880 

0.2120 

.982 

.964 

0.6859 

0.1359 

20-0000* 

1.0878 

0.3522 

.666 

.444 

0.2901 

0.2254 

20-1200 

2.7297 

0.2756 

.957 

.916 

0.4660 

0.1424 

21-0000 

5.0908 

0.3520 

.979 

.958 

0.9154 

0.1960 

21-1200 

5.1735 

0.2328 

.989 

.978 

0.9479 

0.1461 

22-0000 

3.7348 

0.0533 

.999 

.998 

0.6789 

0.0335 

22-1200 

4.5463 

0.1532 

.994 

.988 

0.8320 

0.0965 

23-0000 

6.5325 

0.5259 

.975 

.951 

1.1656 

0.2745 

23-1200 

4.5879 

0.1398 

.996 

.992 

0.8014 

0.0769 

TIME  MEAN 

3.6797 

.892 

STD.  DEV. 

1.4896 

fft(  A  ) 

Corresponds    to   R12     <    -85 


35 


TABLE    IX 

Values   of     X     and    related   statistics 
for   the   data-sample    for  Maniwaki,    Que, 


STANDARD 

CORRELATION 

STANDARD 

STANDARD 

ERROR  OF 

COEFFICIENT 

DEVIATION 

ERROR 

DATE 

X  -VALUE 

D2 

DAY -GMT 

CALCULATED 

X -VALUE 

R12 

R12 

aY 

SYIX 

16-0000 

5.5686 

0.32  28 

.981 

.962 

1.0019 

0.2053 

16-1200 

3.0656 

0.2679 

.971 

.943 

0.4894 

0.1246 

17-0000 

3.6455 

0.1713 

.987 

.974 

0.6561 

0.1097 

17-1200 

3.2185 

0.2048 

.977 

.955 

0.5867 

0.1315 

18-0000 

3.2674 

0.2213 

.976 

.953 

0.6138 

0.1405 

18-1200 

2.8017 

0.1775 

.977 

.955 

0.5597 

0.1249 

19-0000 

4.2382 

0.1084 

.996 

.992 

0.8075 

0.0770 

19-1200* 

1.1753 

0.3235 

.724 

.524 

0.2892 

0.2077 

20-0000 

1.8514 

0.4214 

.874 

.764 

0.3945 

0.2075 

20-1200 

3.0911 

0.1692 

.987 

.974 

0.5533 

0.0945 

21-0000 

3.2247 

0.0641 

.999 

.998 

0.5703 

0.0320 

21-1200 

2.8132 

0.1215 

.991 

.982 

0.5029 

0.0714 

22-0000 

2.3573 

0.2628 

.933 

.870 

0.4481 

0.1680 

22-1200 

2.6480 

0.1180 

.989 

.978 

0.4888 

0.0747 

23-0000 

3.4791 

0.1644 

.989 

.978 

0.6695 

0.1038 

23-1200 

3.8057 

0.3301 

.971 

.943 

0.6607 

0.1670 

TIME  MEAN 

3.1407 

.922 

STD.  DEV. 

0.9849 

«t<  ^  ) 

^Corresponds    to   R,  ~  <    .85 
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TABLE   X 

Values   of     \      and   related    statistics 
for    the   data-sample    for   Sault-Ste-Marie,    Mich 


DATE 
DAY-GMT 

A  -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

A -VALUE 

CORRELATION 
COEFFICIENT 
2 
R12    R12 

STANDARD 
DEVIATION 

STANDARD 
ERROR 

SYIX 

16-0000 

2.4004 

0.3536 

.932 

.869 

0.4486 

0.1741 

16-1200 

2.4229 

0.1205 

.986 

.972 

0.4355 

0.0770 

17-0000 

4.3460 

0.1218 

.995 

.990 

0.7756 

0.0780 

17-1200 

3.6100 

0.1099 

.995 

.990 

0.6639 

0.0697 

18-0000 

3.7061 

0.1273 

.994 

.988 

0.6828 

0.0807 

18-1200 

3.5366 

0.1136 

.994 

.988 

0.6515 

0.0721 

19-0000 

2.0358 

0.2896 

.904 

.817 

0.3889 

0.1733 

19-1200 

2.2482 

0.2289 

.943 

.889 

0.4225 

0.1463 

20-0000 

2.8100 

0.1758 

.979 

.958 

0.5218 

0.1108 

20-1200 

3.2512 

0.0385 

.999 

.998 

0.5917 

0.0242 

21-0000 

3.3835 

0.1394 

.991 

.982 

0.6221 

0.0879 

21-1200 

2.1833 

0.2288 

.945 

.893 

0.4217 

0.1446 

22-0000 

3.0157 

0.1823 

.981 

.962 

0.5614 

0.1153 

22-1200 

3.3318 

0.1228 

.993 

.986 

0.6384 

0.0775 

23-0000 

3.8431 

0.1015 

.997 

.994 

0.7333 

0.0640 

23-1200 

3.7412 

0.1077 

.996 

.992 

0.7496 

0.0680 

TIME  MEAN 

3.1166 

.954 

STD.  DEV. 

0.6940 

M  A  ) 
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TABLE   XI 

Values   of     A      and   related   statistics 
for    the   data-sample    for   Bismarck,    N.    D, 


DATE 
DAY-GMT 

A  -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

A. -VALUE 

CORRELATION 
COEFFICIENT 
2 
R12   R12 

STANDARD 
DEVIATION 

STANDARD 
ERROR 

SYIX 

16-0000 

2.6549 

0.2499 

.955 

.912 

0.5027 

0.1564 

16-1200 

3.4759 

0.2024 

.984 

.968 

0.6213 

0.1180 

17-0000 

6.5097 

0.1724 

.996 

.992 

1.1819 

0.1080 

17-1200 

NO  DATA 

AVAILABLE 

18-0000 

NO  DATA 

AVAILABLE 

18-1200 

3.6856 

0.2127 

.984 

.968 

0.6601 

0.1243 

19-0000 

4.4516 

0.2273 

.986 

.972 

0.8179 

0.1427 

19-1200 

4.5117 

0.1570 

.995 

.990 

0.8724 

0.0955 

20-0000 

3.4355 

0.1859 

.983 

.966 

0.6154 

0.1180 

20-1200 

2.0076 

0.2615 

.918 

.843 

0.3888 

0.1611 

21-0000 

4.9089 

0.1923 

.992 

.984 

0.8086 

0.1088 

21-1200 

4.5735 

0.1341 

.996 

.992 

0.7678 

0.0743 

22-0000 

5.0893 

0.2730 

.987 

.974 

0.9794 

0.1640 

22-1200 

3.8804 

0.0609 

.999 

.998 

0.7341 

0.0382 

23-0000 

3.6326 

0.3285 

.954 

.910 

0.6710 

0.2088 

23-1200 

2.1195 

0.1005 

.987 

.974 

0.3789 

0.0639 

TIME  MEAN 

3.9240 

.960 

STD.  DEV." 

1.2147 

»t<  A  ) 
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TABLE   XII 

Values   of     \      and   related   statistics 
for   the   data-sample   for  Glasgow,   Mont 


DATE 
DAY -GMT 

X  -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

A.  -VALUE 

CORRELATION 
COEFFICIENT 
2 
R12   R12 

STANDARD 
DEVIATION 

STANDARD 
ERROR 

SYIX 

16-0000 

3.7842 

0.1697 

.990 

.980 

0.6406 

0.0943 

16-1200 

6.6178 

0.3483 

.986 

.972 

1.1253 

0.1938 

17-0000 

5.6933 

0.1384 

.997 

.994 

0.9578 

0.0770 

17-1200 

6.0083 

0.2198 

.993 

.986 

1.0154 

0.1224 

18-0000 

6.0189 

0.1978 

.995 

.990 

1.0630 

0.1100 

18-1200 

5.1013 

0.1940 

.993 

.986 

0.8650 

0.1083 

19-0000 

6.0881 

0.1843 

.995 

.990 

1.0284 

0.1028 

19-1200 

2.6928 

0.2074 

.972 

.945 

0.4666 

0.1158. 

20-0000 

3.3309 

0.2818 

.963 

.927 

0.5644 

0.1593 

20-1200 

2.2316 

0.1834 

.965 

.931 

0.3754 

0.1031 

21-0000 

5.3035 

0.2107 

.991 

.982 

0.8700 

0.1187 

21-1200 

3.1120 

0.3134 

.953 

.908 

0.5487 

0.1746 

22-0000 

4.4330 

0.1457 

.995 

.990 

0.7500 

0.0813 

22-1200 

3.3728 

0.1415 

.991 

.982 

0.5573 

0.0802 

23-0000* 

0.6804 

0.3303 

.614 

.377 

0.1954 

0.1649 

23-1200* 

0.8355 

0.2358 

.746 

.557 

0.1885 

0.1316 

TIME  MEAN  4.0815 


906 


STD.    DEV. 


CTt(    A  ) 


1.8749 


^Corresponds    to   R^2     <    -85. 
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TABLE   XIII 

Values   of       A    and   related   statistics 
for   the  data-sample   for   Spokane,    Wash, 


DATE 
DAY -GMT 

\     -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

X  -VALUE 

CORRELATION 
COEFFICIENT 
2 
R12   R12 

STANDARD 
DEVIATION 

aY 

STANDARD 
ERROR 

SYIX 

16-0000 

6.1242 

0.2171 

.994 

.988 

1.0350 

0.1209 

16-1200 

5.3894 

0.0292 

.999 

.998 

0.9060 

0.0163 

17-0000 

6.2015 

0.2579 

.992 

.984 

1.0971 

0.1432 

17-1200 

5.4573 

0.0661 

.999 

.998 

0.9586 

0.0367 

18-0000 

4.1476 

0.2953 

.973 

.947 

0.6969 

0.1673 

18-1200 

3.1086 

0.2754 

.959 

.920 

0.5305 

0.1562 

19-0000 

3.3134 

0.1454 

.990 

.980 

0.5484 

0.0825 

19-1200 

2.5267 

0.1350 

.985 

.970 

0.4199 

0.0766 

20-0000 

2.6138 

0.1618 

.980 

.960 

0.4352 

0.0914 

20-1200 

4.5322 

0.4601 

.948 

.899 

0.7780 

0.2593 

21-0000 

4.2112 

0.3912 

.956 

.914 

0.7166 

0.2204 

21-1200 

3.6604 

0.2030 

.984 

.968 

0.6057 

0.1145 

22-0000* 

1.6059 

0.3356 

.822 

.676 

0.3184 

0.1894 

22-1200* 

1.0266 

0.2085 

.829 

.687 

0.2015 

0.1176 

23-0000 

5.3603 

0.2307 

.990 

.980 

0.8799 

0.1299 

23-1200 

4.0725 

0.4338 

.968 

.937 

0.5157 

0.1406  ' 

TIME  MEAN 

3.9594 

.925 

STD.  DEV. 

1.5414 

M  *  ) 

'^Corresponds    to  Rio  <    «85. 
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TABLE   XIV 

Values   of    A      and    related    statistics 
for    the  data-sample   for  Quillayute,    Wash, 


DATE 
DAY -GMT 

A  -VALUE 
CALCULATED 

STANDARD 
ERROR  OF 

A  -VALUE 

CORRELATION 
COEFFICIENT 
2 
R12    R12 

STANDARD 
DEVIATION 

STANDARD 
ERROR 

SYIX 

16-0000 

5.5314 

0.1401 

.996 

.992 

1.0899 

0.0992 

16-1200 

5.5505 

0.1507 

.996 

.992 

1.1368 

0.1065 

17-0000 

4.3589 

0.2506 

.979 

.958 

0.8494 

0.1789 

17-1200* 

1.4652 

0.3352 

.771 

.594 

0.3626 

0.2394 

18-0000* 

1.7747 

0.3777 

.793 

.629 

0.4275 

0.2701 

18-1200 

1.6247 

0.2586 

.867 

.752 

0.3580 

0.1850 

19-0000 

2.0580 

0.1481 

.968 

.937 

0.4061 

0.1058 

19-1200 

2.3200 

0.1796 

.963 

.927 

0.4596 

0.1282 

20-0000 

2.5319 

0.1793 

.969 

.939 

0.4977 

0.1278 

20-1200 

3.6523 

0.3188 

.954 

.910 

0.7291 

0.2272 

21-0000* 

1.6818 

0.2995 

.842 

.709 

0.3805 

0.2133 

21-1200 

1.4532 

0.3086 

.794 

.630 

0.3485 

0.2199 

22-0000 

1.7184 

0.2594 

.878 

.771 

0.3722 

0.1846 

22-1200 

2.0152 

0.1758 

.982 

.964 

0.4858 

0.1019 

23-0000 

4.1114 

0.1946 

.986 

.972 

0.7917 

0.1382 

23-1200 

3.9581 

0.3402 

.972 

.945 

0.7002 

0.1755 

TIME  MEAN 

2.8629 

.851 

STD.  DEV. 

1.5414 

at(  A  ) 

Corresponds    to   R^2     <     .85 
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In  table  III,  a  low  correlation  of  0.642  occurs  at  22  March,  0000GMT, 
which  results  in  a  relatively  high  standard  error  Sy|x-   Note  in  Table  III 
that  most  of  the  profiles  have  Ri2>0.95  with  associated  small  standard 
errors.   There  were  two  cases  in  Table  III  with  low  correlation  coefficients, 
specifically  0.704  and  0.642.   This  pattern  of  correlation  coefficients 
Ri2  ^  0.85  except  for  a  possible  set  of  one  to  three  profiles  for  which 
R  1 2 f ell  below  this  value,  continues  throughout  all  of  the  12  stations 
considered  here  as  shown  in  Tables  III  through  XIV.   Specific  instances  in 
these  tables  were  the  Ri2~value  value  drops  below  0.85  have  been  marked 
by  an  asterisk  alongside  of  the  date.   In  all  of  these  cases  the  parameter 
X  is  likewise  smaller  than  normal  in  view  of  Eq.  (20) 

X    =  R12(  ay  /  ax)  (44) 

All  of  these  special  cases  where  R]_2  <  0.85  are  characterized  by  a 
standard  deviation  a  y  lower  than  typical  for  the  station.   Such  cases  of 
small   ffY~values  represented  cases  where  statistically  there  was  insufficient 
vertical  variation  in  the  moisture  range  to  ascribe  it  successfully  to  the 
power-profile  parameterization. 

An  example  of  a  smaller- than-normal  correlation  coefficient  R-,2  and 
its  correspondingly  small   X  has  been  extracted  from  the  regression  print 
out  in  the  form  of  Table  XV.   This  table  lists  the  statistical  results 
from  Waycross,  Ga . ,  at  0000GMT,  22  March,  1971.   This  power-profile 
analysis  resulted  in  Ri2  =  0.642  and   X  =  0.9815,  which  was  considerably 
lower  than  the  mean   X  -value.   In  Table  XV,  the  observed  moisture  profile 

is  indicated  in  column  2  as  Y  of  Eq.  (7),  and  in  column  3  is  shown  the 

A 
corresponding  Y  of  Eq.  (40).    The  residual  shown  in  column  4  results  from 

A  A 

forming  Y-Y  from  Eq.  (41).   As  expected,  Y  began  as  a  maximum  positive 

value  at  the  surface  and  decreased  to  0.000  at  500  mb  and  becomes  negative 
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TABLE  XV 

An  example  for  Waycross,  Ga .  (0000  GMT,  22  Mar.,  1971) 
showing  a  condition  leading  to  an  ill-defined    A  -profile 


PRESSURE  (mb) 


Y  (obs) 

A 

Y 

RESIDUAL 

0.4660 

0.3014 

0.1646 

0.3179 

0.2955 

0.0224 

0.2719 

0.2736 

-0.0017 

0.2478 

0.2505 

-0.0028 

0.2020 

0.2262 

-0.0241 

0.1024 

0.2003 

-0.0979 

-0.0060 

0.1728 

-0.1789 

-0.1507 

0.1434 

-0.2942 

-0.1579 

0.1118 

-0.2716 

0.1296 

0.0777 

0.0519 

-0.0082 

0.0406 

-0.0488 

0.0000 

0.0000 

0.0000 

-0.3650 

-0.0449 

-0.3201 

-0.6949 

-0.0951 

-0.5998 

1014 
1000 
950 
900 
850 
800 
750 
700 
650 
600 
550 
500 
450 
400 
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TABLE  XVI 

Environmental  parameters  used  in  the 
calculation  of  Table  XV. 


PRESSURE  (mb) 


RELATIVE 

MIXING 

EMP  (C) 

HUMIDITY,  % 

RATIO  (g/kg) 

16.8 

33 

3.8801 

15.6 

25 

2.7590 

12.0 

27 

2.4816 

8.3 

31 

2.3477 

6.9 

29 

2.1130 

4.3 

26 

1.6800 

1.6 

23 

1.3087 

-   .5 

18 

.9378 

-  2.9 

20 

.9186 

-  6.3 

48 

1.7885 

- 10 . 9 

48 

1.3021 

-16.9 

77 

1.3270 

-19.7 

39 

.5726 

-27.2 

34 

.2678 

1014 
1000 
950 
900 
850 
800 
750 
700 
650 
600 
550 
500 
450 
400 
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above  500  mb .   The  observed  profile,  Y(p),  began  also  as  a  positive  value 
at  the  surface,  but  becomes  negative  in  the  layer  750  to  600  mbs  indicating 
that  the  actual  mixing  ratio  profile  had  values  less  than  W500  in  that  layer. 

Table  XVI  lists  the  actual  atmospheric-sounding  parameters  which  resulted 
in  the  values  shown  in  Table  XV.   Note  that  the  reported  relative  humidities 
are  low  in  the  750  to  650  mb  layer  and  high  in  the  layer  600  to  450  mb. 
Since  by  Eq.  (9),  w  is  a  function  of  relative  humidity  as  well  as  of 
temperature  and  pressure,  relatively  low  w-values exist  from  750  to  650  mb 
and  high  w-values in  the  600-  to  450  mb-  range.   This,  at  least  locally,  is 
inconsistent  with  the  assumption  on  which  the  lower-fit  is  based  and,  as 
calculated,  gives  a  low  value  of  Ri?-   Note  that  the  power-fit  need  not 
have  been  inconsistent  if  the  temperature  lapse  was  large  enough  through 
the  600  -  450  mb  layer,   which  in  this  instance  was  not  the  case. 

As  just  noted,  the  second  major  factor  which  controls  the  consistency 
of  the  profile-fit  at  level  p  is  the  temperature  T(p).   It  can  also  be 
shown,  by  example,  that  a  temperature  inversion  or  an  isothermal  layer 
may  be  (depending  upon  relative  humidity)  associated  with  w  values  which 
do  not  decrease  in  height,  resulting  in  the  low  correlation  coefficients 
R12.   In  Tables  III  to  XIV,  all  cases  where  R-^  -  0.85  which  thus  led  to 
an  ill-defined  (i.e.,  too  small)  value  of   A  are  denoted  by  an  asterisk. 
The  explanation  of  all  of  these  special  cases  is  similar  to  that  advanced 

for  the  Waycross,  Ga .  case  in  Tables  XV,  XVI. 

2 
Note  finally  that  it  is  convenient  to  list  the  time-mean  A.   and  R-^ 

for  each  station  listed  in  Tables  III,...,  XIV.   In  addition  it  is  desir- 
able to  compute  the  temporal  standard  deviation  of  the    A  ' s  in  each 
table  according  to 

V  *  )  = 


(  A.  -  A  )2+...  +  (   A1fi  -  A  )2-^ 


1 


15 


] 


(45) 
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B.   REDUCTION  OF  VARIANCE  OF   A 

Having  obtained  the  sixteen   A  -values  per  station  the  next  portion 
of  the  investigation  involved  determining  the  feasibility  of  reducing  the 
variance  of   A  for  each  station  over  the  period  of  time  considered  using 
other  sounding  parameters.   The  synoptic  variables  introduced  into  the 
computation  as  independent  variables  permitted  utilization  of  the  stepwise 
multiple  regression  aspect  of  the  BIMED  02R  program  where   A  was  the 
dependent  variable.   In  addition  to  the  16  values  of  A   per  station,  16 
values  of  each  of  the  following  X.  were  read  into  the  program  as  indepen- 
dent variables: 

surface  temperature  defined  as  X1 

surface  mixing  ratio  defined  as  X~ 

850  mb  temperature  defined  as  Xo 

850  mb  mixing  ratio  defined  as  X. 

v  4 

700  mb  temperature  defined  as  X,. 

500  mb  temperature  defined  as  X^. 
The  above  variables  were  used  in  the  multiple  regression  as  inputs  with 
the  computed  A   as  the  dependent  variable  in  each  case.   The  desired 
solution  equation  would  then  be  of  the  form  of  Eq.  (32)  where  Y-.  =  A  , 
and  the  X.  are  the  respective  independent  variables.   The  stepwise 
regression  program  BIMED  02R  was  then  utilized  to  derive  the  best-fit 
form  of 

A  =  AQ  +  ALXL  +  A2X2  +. . -+A6X6  (46) 

Among  the  factors  which  might  be  considered  as  useful  in  reducing  the 
temporal  variance  of  A  are  variations  in  the  synoptic  patterns,  and  those 
which  may  reflect  purely  local  conditions.   Among  the  latter  effects 
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(e.g.,  Tables  III,...,   XIV),  station  location  and  elevation  should  be 

considered,  as  well  as  proximity  to  large  water  bodies  and  lakes.   In 

the  sounding  variables  X- , . .  .  ,X  ,  the  lower-level  parameters,  were  con- 

*■  6 

sidered  to  be  primarily  representative  of  local  effects,  whereas  the 
higher-level  variables  were  presumed  to  be  indicative  of  the  larger-scale 
synoptic  time-scales,  although  there  is  undoubtedly  some  interaction 
between  the  low-  and  high-level  variables. 

The  temporal  standard  deviations  o    (   A  )  listed  at  the  base  of 


of  <7t(  A  )  =  0.694  at  Sault  Ste  Marie,  and    O     (  A.   )  =  0.701 


Tables  III,..., XIV  resulting  from  Eq.  (45)  ranged  from  minimal  values 

at 

San  Diego  to  a  maximum  (J    (      A   )  =  1.875  at  Glasgow,  Montana.   In  general 
it  can  be  seen  that  stations  close  to  large  bodies  of  water  (Charleston, 
San  Diego,  Maniwaki,  and  Sault  Ste  Marie)  have  relatively  small  standard 
deviations  whereas  stations  away  from  bodies  of  water  have  larger  standard 
deviations  (Bismarck,  Glasgow,  and  Spokane).   The  cases  of  large  values 
of  a  .(   A  )  seem  to  afford  considerable  opportunity  for  utilization 
of  the  regression  procedure  for  the  reduction  of  the    A  -variance.   How- 
ever, it  was  decided  to  test  the  diagnostic  usefulness  of  Eq.  (46)  for  all 
twelve  stations,  regardless  of  the  o     (  A   )  value. 

A  preliminary  test  run  was  conducted  using  the  BIMED  02R  program  with 
the  variables  defined  previously  in  the  section.   During  this  run  all 
independent  variables  were  allowed  to  enter  freely  based  solely  on  their 
partial  correlation  values.   The  program  was  allowed  to  run  until  the 
requirement  of  2%  added  fractional  explained  variance  after  introduction 
of  the  kth  variable  could  not  be  met. 

Two  facts  became  evident  from  this  preliminary  test:   (1)  the  indepen- 
dent variable  with  the  highest  initial  correlation  coefficient  in  almost 
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all  cases  was  X^,  and  (2)  the  maximum  number  of  significant  steps  necessary 
to  ensure  the  required  minimum  added  explained  variance  at  each  step  was 
four.   At  steps  five  and  six  of  the  preliminary  run  the  added  explained 
variance  upon  predictor-entry  was  too  small.   Since  a  low-level  moisture 
predictor  was  desired  to  be  compatible  with  the  variations  of  the   X  - 
profile  parameter  (cf.,  the  discussion  for  Waycross  at  0000  GMT,  22  March, 
1971),  it  was  decided  to  utilize  the  control-delete  option  of  BIMED  02R 
to  force  in  Xa  as  the  first  predictor  at  each  station.   The  second  option 
used  was  that  of  limiting  the  number  of  predictors  to  four  in  each 
regression  analysis.   The  results  of  the  sets  of  regressions  over  the  12 
stations  are  shown  in  Table  XVII. 

Equation  (42)  gives  the  relationship  between  the  standard  error  of 
X  ,  after  the  four-step  multiple  correlation  R  has  been  determined 

S  AjX  =  (1  -  h£)%   V  A  )  (47) 

From  Table  XVII  it  can  be  seen  that  when  a    (  A  )  is  relatively  small, 
the  explained  fraction  of  this  statistic  is  relatively  small.   The  result, 
then,  of  application  of  the  multiple  regression  technique  is  to  reduce 
the  standard  deviation  Q    (  A   )  of  all  stations,  but  to  reduce  the 
stations  with  the  largest  O    (  \    )  the  most. 

San  Diego  is  a  good  example  of  a  station  having  a  small  o    (     \    ), 
°    (  A  )  =  0.7012,  which  means,  according  to  the  above  result,  a 
relatively  small  reduction,  specifically  to  S\iy=0.5686.   Correspondingly, 
this  small  additional  reduction  is  reflected  in  the  relatively  low  multiple 

correlation  coefficient  at  San  Diego  of  R  =  .5852.   At  the  other  extreme 

°  m 

is  Glasgow  having  a  temporal  standard  deviation  O     (  A  )  =  1.8749.   By 
the  multiple  regression  procedure  a  standard  error  of  estimate  of  0.8745 
results.   Again  this  is  consistent  with  the  high  multiple  correlation 
coefficient,  E^,  =  0.8846. 
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Figure  2:   Histogram  depicting  number  of  stations  within  certain 
standard  error  limits. 


50 


The  result  of  the  multiple  regression  for  each  station  is  to  give  a  re- 
duced standard  error  of  A   .   A  histogram  depicting  the  12  standard  error 
values  is  shown  in  Figure  2,  on  a  scale  of  S  \ i  x  from  0.50  to  0.95.   Note 
that  most  of  the  standard  errors  after  regression  fall  into  the  range 
0.5  to  0.65. 

The  righthand  column  in  Table  XVII  lists  the  order  of  entry  of  the 
independent  variables  into  the  stepwise  regression  and  is  based  on  the 
value  of  the  partial  correlation  of  the  respective  variables  as  explained 
in  Section  III  (B) .   From  Table  XVII  it  may  be  noted  that  in  seven  of  the 
12  stations  the  second  independent  variable  added  to  the  regression  was  X,-. 
Also  in  seven  of  the  12  stations  the  third  independent  variable  added  was 
X2  •   This  led  to  the  idea  of  possible  classification  of  stations  according 
to  a  forced  order  of  entry,  an   option  of  BIMED  02R,  and  will  be  discussed 
in  the  next  sub-section. 

C.   STATION  CLASSIFICATION  ACCORDING  TO  THE  PREDICTOR  ENTRY 

Table  XVIII  lists  an  attempt  to  classify  stations  by  the  relative 
priority  of  the  independent  variables  to  the  program  specified  in  the  order 
X4,  X/-,  X  ,  and  the  remaining  predictor  in  a  4-step  selection  left  as  free 
choice  according  to  the  selection  criterion.   As  can  be  seen  by  comparing 
Tables  XVII  and  XVIII  there  was  no  change  in  the  multiple  correlation 
coefficient  at  10  of  the  12  stations  as  a  result  of  specifying  the  4-6-2 
input  order.   This  is  because  in  almost  all  cases  the  variables  X^,  X^,  and 
X2  were  in  the  original  run  (Table  XVII)  even  though  not  necessarily  in  this 
specific  order.   Table  XVIII  lists  the  "single  free  choice"  variable  under 
order  of  input  as  the  last  independent  variable  in  the  list  for  each  station 
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It  may  be  noted  that  there  are  four  stations  with  surface  temperature  as 
the  last  independent  variable  (Charleston,  Maniwaki,  Bismarck,  and 
Spokane),  four  with  500  mb- tempera ture  (Waycross,  Lake  Charles   El  Paso, 
and  Caribou),  and  four  with  850  mb-temperature  (San  Diego,  Sault  Ste 
Marie,  Glasgow,  And  Quillayute) .   There  might,  then,  be  some  way  of 
classifying  or  typing  stations  on  the  basis  of  the  free-choice  predictor 
selected  at  step  4.   However,  with  the  small  station  set  available   it 
was  not  clear  that  this  could  be  done. 

Finally,  it  is  to  be  emphasized  that  Table  XVII  indicates  that  the 
observed    A.  -values  tend  to  adjust  to  (a)  the  synoptic-scale  time 
variations  in  the  sounding  and  to  (b)  the  local-effect  variations.  The 
latter  effect  is  harder  to  specify  so  that  one  achieves  some  limiting 
accuracy  in  the  expected  value  resulting  from  the  regression.   The  magnitude 
of  this  accuracy  is  suggested  by  the  standard  error  which  tends  to  show  a 
cluster  of  values  in  the  range  0.50  to  0.65  (Fig.  2). 
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V.   PRECIPITABLE  WATER  VAPOR  TESTS 

As  noted  in  the  introduction,  two  verification  techniques  were  to  be 
used  to  evaluate  the  results  of  this  investigation.   The  first  of  these, 
namely  verification  of  the  moisture  profile,  has  been  discussed  in  detail 
in  Section  IV  from  a  statistical  viewpoint.   In  this  section,  an  analysis 
of  the  verification  of  the  moisture-parameterization  technique  from  the 
viewpoint  of  precipitable  water  vapor  is  presented. 

The  formulas  used  to  calculate  the  precipitable  water  are  shown  as 
Eqs.  (12)  and  (13).   Equation  12  was  utilized  to  obtain  precipitable  water 
based  upon  the  moisture  profile  of  the  data  for  the  given  data  sample.   It 
is  the  precipitable  water  from  the  observed  sounding,  and  hence  will  be 
referred  to  as  observed  precipitable  water,  WQ.   Equation  13  is  used  to 
calculate  the  precipitable  water  utilizing  the    A  -values  already  obtained. 
As  can  be  seen  from  Eq.  (13),  it  is  based  on  only  three  variables:   500-mb 
mixing-ratio,  surface  pressure,  and  X   .   There  are,  however,  two  possible 
methods  since   A  -values  were  obtained  by  the  profile-parameter  method 
and,  secondly,  by  the  multiple  regression  technique.   Precipitable  water 
calculated  using  the  profile  parameter   A  -values,  one  for  each  sounding, 
will  be  designated  as  Wcs  and  those  relatively  smooth  values  obtained  by 
using  the  multiple  regression   A  -values  denoted  by  W   . 

For  this  phase  of  the  investigation,  four  stations  were  chosen  primarily 
on  the  basis  of  their  varying  geographical  locations:   Waycross,  Ga . ,  in 
the  southeastern  United  States;  San  Diego  in  the  southwest;  Spokane,  Wash., 
in  the  northwest;  and  Sault  Ste  Marie,  Mich.,  in  the  northeast.   Results 
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of  the  precipitable  water  calculations  are  listed  in  Tables  XIX  through 
XXII  and  are  plotted  in  Figures  3  through  6. 

As  can  be  seen  from  Tables  XIX, ..., XXII  or  from  Figures  3,..., 6,  the 
W_c-value  is  in  almost  all  cases  within  approximately  0.05  cms  of  the  W. 
value.   Needless  to  say,  with  such  good  verification,  the  trend  of  W 
is  the  same  as  the  observed  trend.   It  is  obvious  from  the  results,  then, 
that  the  use  of  the  power-profile  technique  can  be  expected  to  yield 
precipitable  water  values  which  are  very  representative  of  the  WQ  values, 
i.e.  within  0.05  cms  (cf.,  Tables  XIX, . . . ,XXII) . 

The  values  of  Wcm  are  also  presented  in  Tables  XIX, ..., XXII  and  in 
Figures  3,..., 6.   As  can  be  seen  errors  do  exist  in  the  Wcm  values, 
assuming  WQ  is  the  correct  value.   These  errors  can  at  times  become  large, 
i.e.  as  much  as  100%.   In  general,  however,  the  Wcm  values  which  result 
from  using   A  -values  obtained  by  multiple  regression  follow  the  time 
trend  of  W  .   Large  errors  in  Wcm  are  then  the  exception  rather  than  the 
rule,  and  as  can  be  seen  in  Figures  3,..., 6,  occur  at  only  one  or  two 
sounding  times  per  station.   The  soundings  which  give  the  largest  Wcm 
errors  are,  in  fact,  just  those  soundings  which  had  the  lower  values  of 

A  and  of  the  correlation  coefficients  R^2  i-n  connection  with  the  Wcs 
computations . 

The  0000  GMT,  22  March,  1971  sounding  at  Waycross,  Ga .  will  be  used 
to  illustrate  the  point  of  a  large  WCm  error.   As  can  be  seen  in  Table  III, 
this  sounding  gave  a  low  correlation  coefficient  value,  0.642,  and  a   A  - 
value  of  0.9815  where  the  time-mean  over  the  16  soundings  for  that  station 
was  2.8557.   Figure  3  shows  that  the  use  of  the    A  -value  from  the  power- 
fit  profile,   A  =  0.9815,  in  Eq.  (13)  gives  an  accurate  verification  of 
W0.   On  the  other  hand  the  value  of   A  resulting  from  the  multiple 
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TABLE  XIX 


Table  of  precipitable  water  vapor  values  (W  ,  W  s,  Wcm) 


for  Waycross,  Ga. 


DAY 

OBSERVED 

PROFILE 

MULT.  REG. 

DAY  -  GMT 

Wo 

W 

.  cm 

16  -  0000 

3.7190 

3.8704 

3.0169 

16  -  1200 

1.0227 

1.1884 

1.2003 

17  -  0000 

0.9360 

0.8618 

0.4918 

17  -  1200 

0.5325 

0.5421 

0.5940 

18  -  0000 

0.5678 

0.5863 

0.5923 

18  -  1200 

0.6549 

0.6600 

0.5514 

19  -  0000 

1.1639 

1.1672 

1.6403 

19  -  1200 

2.2758 

2.2377 

1.9609 

20  -  0000 

1.0744 

0.9818 

0.6177 

20  -  1200 

0.5083 

0.5216 

0.5058 

21  -  0000 

0.4897 

0.4766 

0.5948 

21  -  1200 

0.8447 

0.8370 

0.9129 

22  -  0000 

0.9807 

1.1674 

2.4002 

22  -  1200 

1.2914 

1.1747 

1.1335 

23  -  0000 

1.5845 

1.3996 

0.9415 

23  -  1200 

3.0389 

2.8514 

2.1848 
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TABLE  XX 

Table  of  precipitable  water  vapor  values  (W0,  Wcs,  Wcm) 
for  San  Diego,  Calif. 


DAY 

OBSERVED 

PROFILE 

MULT.  REG. 

DAY  -  GMT 

Wo 

WCs 

Wcm 

16  -  0000 

0.7488 

0.7849 

0.5459 

16  -  1200 

0.8116 

0.8317 

0.6233 

17  -  0000 

0.8222 

0.8833 

0.6292 

17  -  1200 

0.8586 

0.8796 

0.8950 

18  -  0000 

0.8783 

0.8523 

0.7305 

18  -  1200 

1.0573 

1.0334 

0.7619 

19  -  0000 

0.9586 

1.0137 

0.8904 

19  -  1200 

1.1637 

1.2238 

1.6504 

20  -  0000 

0.7991 

0.9855 

1.5071 

20  -  1200 

0.7310 

0.8049 

0.8771 

21  -  0000 

0.6700 

0.7594 

0.8141 

21  -  1200 

0.9416 

0.9233 

0.8065 

22  -  0000 

0.7205 

0.8268 

1.0792 

22  -  1200 

1.2398 

1.3404 

1.8782 

23  -  0000 

1.1738 

1.3048 

1.7369 

23  -  1200 

1.7495 

1.7309 

1.4738 
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TABLE  XXI 

Table  of  precipitable  water  vapor  values  (W0)  Wcs,  Wcm) 

for  Spokane,  Wash. 


DAY 

OBSERVED 

PROFILE 

MULT.  REG. 

DAY  -  GMT 

W0 

Wcs 

wcm 

16  -  0000 

0.3970 

0.3932 

0.2822 

16  -  1200 

0.3309 

0.3313 

0.2976 

17  -  0000 

0.3965 

0.4002 

0.4135 

17  -  1200 

0.3807 

0.3731 

0.6097 

18  -  0000 

0.3491 

0.3174 

0.2582 

18  -  1200 

0.4091 

0.3976 

0.3855 

19  -  0000 

0.4692 

0.4820 

0.4334 

19  -  1200 

0.5359 

0.5436 

0.5835 

20  -  0000 

0.5297 

0.5337 

0.6059 

20  -  1200 

0.6018 

0.4964 

0.3805 

21  -  0000 

0.4113 

0.3489 

0.3469 

21  -  1200 

0.3094 

0.2940 

0.1853 

22  -  0000 

0.3098 

0.3036 

0.4224 

22  -  1200 

0.7075 

0.7759 

0.9983 

22  -  0000 

0.7988 

0.7595 

0.6473 

23  -  1200 

0.9039 

1.0635 

1.3750 
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TABLE  XXI 

Table  of  precipitable  water  vapor  values  (W0,  Wcs,  Wcm) 
for  Sault  Ste  Marie,  Mich. 


DAY 

OBSERVED 

PROFILE 

MULT.  REG. 

DAY  -  GMT 

W0 

WCs 

wcm 

16  -  0000 

1.1656 

1.1476 

1.0986 

16  -  1200 

0.3112 

0.3207 

0.3864 

17  -  0000 

0.3320 

0.3293 

0.1813 

17  -  1200 

0.3511 

0.3401 

0.3188 

18  -  0000 

0.3865 

0.3693 

0.3514 

18  -  1200 

0.4089 

0.3887 

0.2807 

19  -  0000 

0.5890 

0.5605 

0.7646 

19  -  1200 

0.7004 

0.7405 

0.8876 

20  -  0000 

0.7978 

0.7349 

0.6034 

20  -  1200 

0.5112 

0.5162 

0.4981 

21  -  0000 

0.5046 

0.5150 

0.4953 

21  -  1200 

0.2910 

0.2889 

0.4080 

22  -  0000 

0.3454 

0.3433 

0.4346 

22  -  1200 

0.2752 

0.2701 

0.2942 

23  -  0000 

0.3232 

0.3234 

0.3215 

23  -  1200 

0.2074 

0.2032 

0.1921 
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regression  procedure  for  the  date  and  station  considered  was  X    -   2.6793. 
Since  the  latter  value  is  considerably  larger  than  the  sounding  best-fit 
value  of   A  ,  the  application  of  the  multiple  regression  \     -value  in 
Eq.  (13)  gives  a  much  larger  value  of  Wcm  than  W  ;  hence  the  large  error, 
W   -  W0,  at  this  sounding  time  (Figure  3). 

The  results  of  Tables  XIX,..., XXII  (and  the  four  supporting  diagrams 
Figures  3,..., 6)  have  been  presented  in  the  statistically  summarized  form 
as  Table  XXIII.   In  Table  XXIII,  the  mean  W  -values  have  been  compared 
to  the  computed  mean  values  of  w*cs  and  W   using  the  16  sounding  samples 
at  the  four  stations,  Waycross,  San  Diego,  Spokane,  Sault  Ste  Marie. 
The  results  were  by-products  of  the  correlation  of  W   separately  against 
W   ,  and  of  w*0  with  W"cm,  which  was  applied  to  the  two  sets  of  data  in 
Tables  XIX, . . . ,XXII,  respectively. 

The  standard  deviations  of  W0,  Wcs,  and  Wcm  have  been  summarized  in 
Table  XXIII  (column  3).  Note  that  the  means  and  standard  deviations  of 
the  W-samples  computed  by  the  different  procedures  are  nearly  equal. 

The  correlation  coefficients  R(WQ,  Wcs)  and  R(WQ,  Wcm)  are  listed  for 

the  four  stations  in  column  4  of  Table  XXIII.   The  corresponding 

2 
coefficients  of  determination  R  are  listed  in  column  5,  and  are  to  be 

interpreted  as  the  fractional  explained  variance  (of  W  ) .   Finally  the 

column  "Std.  Error"  is  to  be  interpreted  as  the  mean  dispersion  of  the 

curve  of  W_  minus  W-c,  and  also  that  of  W_  minus  W„m,  using  the  following 

equation  for  standard  error: 

Swo|Wc  -  (1  -  Re) k      °MQ  (48) 

Note  that  the  estimation  of  Wcs  to  give  W   leads  to  a  standard  error  (a 
mean  residual)  of  0.07  cm,  whereas  the  corresponding  standard  error  in  the 
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comparison  of  Wcs  with  W0  is  close  to  0.31  cms.   This  larger  error  arises 
primarily  from  the  southerly  stations  in  Table  XXIII. 

In  conclusion,  the  value  of  Wcs  based  upon  a  best-fit   X  -value  in 
conjunction  with  Eq.  (13)  gives  a  good  estimation  of  W  .  On  the  other  hand, 
the  computation  of  Wcm  wherein  the   A.  -value  is  estimated  by  a  multiple 
regression  formula  using  gross  parameters  from  the  sounding  has  some  skill 
in  estimation  of  W  ,  but  perhaps  needs  additional  study  for  selection  of 
input  variables  for  optimal  inclusion  in  the  multiple  regression  equation 
used  to  determine   A  . 
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VI.   CONCLUSIONS 


At  a  group  of  12  geographically  diverse  stations,  the  mixing-ratio 
profile  was  parameterized  as  a  power-law  profile,  Eq.  (6),  having  the 
characteristic  exponent   A  .   For  each  sounding   A  was  derived  by  a 
least-squares  fit  with  an  accuracy  usually  characterized  by  a  correlation 
coefficient  of  0.95  or  higher.   At  each  station,  the  fit  was  most  accurate 
in  this  form  when  the  standard  deviation  of  observed  mixing-ratio  was 
large.   The  set  of   A  -values  were  grouped  together  in  time  series  at  each 
station,  and  they  were  found  to  be  quite  variable  in  general.   This 
variability  of  A  was  related  to  the  detailed  structure  of  the  moisture 
profile  from  the  sounding  data.   A  sizeable  part  of  the  standard  devia- 
tion <r  (  A  )  was  found  to  be  explainable  by  use  of  a  multiple  regression 
procedure  which  incorporated  easily  defined  properties  of  the  sounding. 
The  primary  test  of  the  A  -profile  parameterization  was  made  in  application 
to  the  comparison  of  observed  versus  computed  values  of  the  precipitable 
water  vapor  at  four  stations,  with  the  latter  given  by 

(W500)(500) 


w 


[CO   f/Ps    \     A  +1         /Ptop\  A  +1 1 

i)  L\500/  \  500/ 


g(  A  +  1)  IV  500/         \  500  /  (13) 

On  an  individual  sounding  basis  the   A  -value  gave  good  verification 
of  the  observed  precipitable  water  vapor  value.   Over  a  period  of  time 
it  was  possible  to  generate  multiple  regression  equations  for  A  which 
gave  moderately  good  verification  when  applied  for  the  appropriate  sound- 
ing using  Eq.  (13).   There  is  some  possibility  that  this  procedure  can  be 
employed  in  a  predictive  sense  provided  the  parameters  used  in  the 
regression  can  be  successfully  generated  by  a  primitive  equations  model. 
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It  is  possible  that  the  multiple  regression  sounding  parameters  may  be 
selected  to  yield  closer  agreement  with  the   \  -values  determined  from 
the  individual  soundings. 
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